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ABSTRACT 

Iterative Extended ^ickel is modified by inclusion of explicit effective 
internuclear and electronic interactions. The one-electron energies 
are shovm to obey a variational principle because of the form of the 
effective electronic interactions. The modifications permit mimicking 
aspects of valence bond theory with the additional feature that the 
energies associated with valence bond-type structures are explicitly 
calculated. In turn, a new hybrid molecular orbital-valence bond 
scheme is introduced which incorporates variant total molecular elec- 
tronic density distributions similar to the way that Iterative Extended 
Hiickel incorporates atoms. Nitrobenzene is used to illustrate the 
approach. 
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1. Introduction 


Iterative Extended Hlickel (hereafter lEH) has been successfully applied 
to complex and diverse systems. These systems have ranged from por- 
phyrins [1,2] to silica model sites [3]. lEH Is one of the few semlem- 
plrlcal appr.iaches which can Inexpensively deal with large systems 
(e.g., greater than 40 atoms, 150 electrons). Moreover, most of the 
first 92 elements can be readily Incorporated Into an lEH program. As 
is true with any semiempirical method, there are inherent deficiencies 
in lEH [4,5]. Only grocss comments can be made about geometries. For 
example, lEH does not explicitly treat nuclear-nuclear Interaction 
terms. In fact, It does not incorporate either an effective nuclear- 
nuclear algorithm or an electronic-nuclear algorithm. lEH is an "atoms- 
in-molecule" approach and achieves a measure of sensitivity to the 
molecular environment through the charge self-consistency requirement. 
The self-consistent charge requirement does not yield a lower-bound 
variational energy for lEH, however. That is, the total electronic 
energy, which is a negative quantity, decreases in magnitude as charge 
self-consistency is achieved. Thus, an energy price is exacted in the 
lEH scheme by delocalizing electronic density. Useful modifications 
would correct the deficiencies mentioned above while retaining lEH's 
speed and versatility. It is the purpose of this note to report that 
such modifications (which also enhance sensitivity to the molecular 
environment) indeed are possible and can be incorporated with minimal 
disturbance to the desirable features of lEH. These modifications are 
made within the framework of the lEH "philosophy." That is, the struc- 
ture of lEH is semiempirical, and all modifications are made within 
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that structure. The end product appears to be a much more versatile 
semiempirical program. It avoids or compensates for the pitfalls asso 
dated with standard lEH or GNDO, while successfully dealing with sys- 
tems that lEH or CNDO cannot cope with. (This is detailed in the 
accompanying report.) It is still a semiempirical program, with con- 
comitant limitations, that will be described in this note. We begin 
with a brief review of the relevant features of lEH. 


2. Review of Standard lEH 

lEH offers a prescription for the Hamiltonian matrix elements. The 
diagonal elements of the Hamiltonian matrix for the A + 1 iteration 
are given by 


^k(i)ik(i) “ “ l‘^i,,^l^^k(i) l'*i,JlI\(i)' 


( 2 . 1 ) 


Here, q, „ is the charge on the 1th atom from the ith iteration; 
n 

E^(i), *^he valence orbital ionization energies of the kth 

orbital on the i^th atomic center for the neutral or singly 

charged center. The proper is determined by q^, If 

q^ ^ < 0, then chosen, for example. The prescription for the 

off-diagonal elements does vary [6], but the Wolfberg-Helmholz formula 
is representative, where 


,(ii+l) 

k(l),n(j) 
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(m) 

k(i) ,k(i) 


+ H 


(i+1) 

n(j),n(j) 



®k(i) ,n(j) ■ 


( 2 . 2 ) 


Here, K is an interaction constant, usually taken to be greater than 
1.75. In porphyrin work, it is set equal to 1.89 [1], for example. The 
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^k(ll^k(i)' *^n(j)^n(j) diagonal matrix elements for the kth 

orbital on the Uh atom, and the nth orbital on the j^ch atom in the 

£ + 1 iteration. The are the integrated overlap between 

the k(i) and n.(j) orbitals. The wavef unctions used to determine 

S, are Slater-type orbitals optimized to agree with SCF over- 

k(i),n(j) 

lap values, at least for bonding regions. Charge density is calculated 
using the Mulliken approximation where, after the £th iteration, the 
charge 


,<m) 


^i.Jl “ ^1 


(valence) ‘ (E | E ''t <=t!k(i) 

\ h L J 


k(l) 

n s 1) 

\^t,k(i)®t,r(j)^k(i),t(j)|j. 


(2.3) 


+ E 

k(i) ,r(j) 
ii^j 

( 0 ) 

The coefficients k(i)* example, are the coefficients appearing 

in the t molecular orbital wavef unction, and they are associated with 

the k(l) orbitals in the ith atom. Again, S, are the inte- 
- ~ k(i),r(j) 

grated overlaps between the orbitals on the ^ and ^ atomic centers. 

The Hj. are the occupation numbers of the molecular orbitals. The term 


(valence) 


Z. 

X 


p (valence) -1 


p=l 


(2.4 


where Z^ is the nuclear charge on the ^th atom and the sum is over all 
shells below the valence shell. The one-electron molecular energies 
are assumed to contribute to the total energy. 


E 


total 


E ”t‘ 


(2.5) 
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3. The Modlgicatlon 


It has been argued using Extended HUckel (and by extension lEH) that 
nuclear-nuclear interactions are implicitly contained in to 

some degree [4]. Certainly the degree is indeterminate. The state of 
affairs is amended by assuming that the electronic density found after 
the _£th iteration plays a role in determining explicit effective 
nuclear-nuclear intarection terms. The effective internuclear inter- 
action between any two atomic centers is calculated by introducing the 
artifice that each center is surrounded by a non-conducting sphere whose 
radius equals one-half the distance between the centers. The Mulliken 
approximation is adopted along with the approximation that the probable 
electronic density within a sphere is determined by the radial portion 
of the atomic wavef unctions . The latter approximation insures Invari- 
anca o£ the results to special orientation. Thus, if is the 

distance between the j^th and jLth atoms, the product of the screened 
nuclear charges (for the two atomic centers), after the ^th iteration, 
is given by 


(Z Z ) 

i^j ''effective 


= (Zi(valence) - 

^ £ \^S(i)^S(j)^k(i),r(j)|^ 


,k(x) 


k(i) ,r(j) 




(valence) - ^t!i(j) 


r(j) 


Is 0 


+ h 


L \^t,k(i)^t,r(j)\(i),r(j)|) • 
k(i),r(j) 

±l^j 


1 ) 


The and 

for the k(i) and _r(j) 

two centers Increases, 


are the radial portion of the STO functtions adopted 
atomic orbitals, As the distance between the 
the term 


/'/ 7 \ (^) 

'•V'j’' effective 


->• 





op , 


(3.2a) 


That is, the ej.nression in Eq. (3.1) approaches the product of the net 
charges on the i.th and j^th atoms after the ^th Iteration at large dis- 
tances. On the other hand, as the atoms are brought together, 


<2l^j>Sectlve Zi(valence) (valence) , R^. 0. (3.2b) 

The effective nuclear-nuclear interaction term after the _£th Iteration 
is given by 

g < 3 . 3 ) 

in atomic units. Consequently, 


E<^) 

total 


= En,c<« 


+ J(ZiZj) 


a) /R 
effective' xj ’ 


(3.4) 


The second type of modification attempts to address the problem of how 
an electron in an atomic orbital on a particular atomic center "experi- 
ences" its molecular environment. The valence orbital ionization 
energies (hereafter, VOIE) implicitly do account for one-center two- 
electron interactions to a large degree. Requiring charge self- 
consistency does result in the molecular environment impinging on an 
atomic center and, in turn, being affected by that center. However, 
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standard lEH does not explicitly treat multicenter electronic interac- 
tions. Such interactions are explicitly included in the raodiflcatlon 
scheme in an approximate fashion. Since the extent to which the VOIE 
and self-consistent charge requirement compensate for one-center and 
multicenter electronic interaction terms is uncertain, the modification 
will be sufficiently flexible that this very uncertainty can be used as 
a minimization condition. Consider an electron localized in an orbital 
on the _lth atomic center, and assume the validity of the Mulliken approx- 
imation; an expression for the coulomblc Interaction between this elec- 
tron and the electron density on the j^th center then can be set equal to 


fie 


(Z) 

k(i),j 


i*j 


' S Pi(J)/‘*’k(i)''k(i),n(d)^n(j)'^''* 
n(j) 


(3.5) 


Here, a. . is a parameter dependent only upon the order and nature of 

J 

the pair of elements involved in expression in Eq. (3.5), e.g., carbon, 
hydrogen. As implied, it is not assumed that a . . = a, . . The elec- 

J ,1 

( 0 ) 

tronlc density refers to the electronic density assigned to the 

^th orbital on the j,th atomic center after the J^th iteration. All one- 
center terms are suppressed, since the VOIE are chosen to compensate 
for one-center terms. The total coulombic-type Interaction of a 
localized electron in the kth orbital, _lth atomic center in a molecular 
environment after the _£th iteration, with reference to Eq. (3.5), is 
given by 


fie 


(Z) _ 

k(i) “ 


E fie 
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(Z) 

k(l) ,j 


(3.6) 
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In order Co ensure rotational invariance, the <j)B appearing in the 

integral in Eq. (3,5) are restricted to the radial portion of the STO. 

That is, all coulombic interactions are treated as though occurring in 

(i) 

a spherically symmetric environment. The modifies the 

Hamiltonian matrix elements which Involve the kth orbital on the j^th 
center. With reference to Eq. (2.1), the element 


H, 


a+1) 

k(i),k(i) 




i.i' k(i) 


+ <Se 


w 

k(l)’ 


(3.7) 


( 0 ) ( 0 ) 

Consequently, the off-diagonal elements “ n(j) k(l) 

will be modified (following the prescription in Eq. (2.2)). It will be 
shown that, if a. . > 0 for all atomic types, the resulting set of 
molecular energies obey a variational principle. 


In standard lEH, < 0 


,(m) 

'k(i),k(i) 


k(i) 


’]■ 


= 0 , in general. The diagonal 


elements of the Hamiltonian matrix with “ 0 will be denoted by 

^k(i) \(i) j ^ '^^k(i) “ ® ^ since the 

coulombic integrals are positive quantities as are the electronic den- 
sities. The equality holds only if the _ith atom is infinitely distant 
from all other atomic centers. Therefore, for all cases of interest 
'^^k(l) ^ term decreases the absolute magnitude of ^j^(i)\(i) 

with respect to that of turn, the absolute magnitude 

of ^k(i)^n(j) decreased with respect to ^^(i)^n(j) 

assumption 


<a)!k(i)«) ^ 0 


and the additional constraint that 


0 < 


( 1 ) 

k(l) 


< 111 


a+1) 

k(l) ,k(l) 


( 0)1 


for any k(l) , enables us to claim that when the mth eigenvalue of the 
Hamiltonian matrix < 0, Chen 



(3. 8a) 


upon appeal to variational arguments (6, 7] 
i(f+l)(0) > 0, 


Stollarly, when 


m 




, (Ji+1) 
"m 


(3.Sb) 


from variational arguments [7; 2, Theorem 17] . Thus, viewing the 
Hamiltonian matrix as being parameterized by the set of whose 

elements are positive and whose effect Is to decrease the magnitude of 
the elements of the Hamiltonian matrix, yields the inequalities in 
Eqs. (3.8a,b). A brief overview concerning the Interrelationship be- 
tween the VOIE and the electronic interactions (basically coulomblc in 
this scheme) might prove useful. The decrease the VOIE (which 

Is assumed to be a negative quantity). As a consequence, an electron 
will be held less tightly in that orbital on the _lth atomic center. 
Formally, a smaller VOIE for the k(i) orbital can be shown to result 
in a decreased participation of k(l) in any particular bonding MO 
than would be the case if were zero. In turn, the electronic 

density assigned to k(i) by the Mulliken approximation will be smaller 
and, hence, the charge on 1 larger. The electron density that must 
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be rsdlstrlbiited in order to conserve the total number of electrons 
will be redistributed among the remaining orbitals In the system. The 
additional electronic density will decrease the magnitudes of the VOIE 
of these orbitals, with reference to Eq. (2.1) and the following 
inequalities! 


E. 


Ic(i) %(1) " ‘^k(l) 


0 


< E,. 


pO 

*'k(i) 


^ 0 1 


These inequalities are obeyed for any jc(i). Consequently, the total 
molecular energy of the system is decreased in magnitude. 


4. An Example of Imposed GharRe Separation 

Nitrobenzene was chosen as a test case for several reasons! (1) It is 
sufficiently complex that raultlcenter processes can be readily examined, 
yet it is an uncomplicated aromatic system. This simplicity is reflected 
by the fact that naive valence bond arguments correctly predict the 
attack site for electrophilic reaction, as an example. (2) It was of 
Interest to explore what a blend of simple valence bond concepts with 
lEH would yield. All calculations were performed on the NASA-Aines 
CDC 7600. Firstly, the nitrogen electronic interaction terms were 
adjusted so that a stable charge of ca. +1 on the N was obtained when 
charge self-consistency was achieved. No other atomic center had 
explicit electronic interaction terms. This case was labelled NITRO in 
Fig. 1. The next case considered was inclusion of electronic interac- 
tion terms on the ortho carbons while maintaining a charge of ca. +1 
on the nitrogen. It was found necessary to include both ortho carbons; 
asymmetric inclusion of electronic Interaction terms resulted in 
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oBolllatlon. Similarly, charge self-conalstency could not be achieved 
for charges much greater than +% on each of the ortho carbons. This 
case is labelled ORTHO in Fig. 1. Finally, the para position carbon 
was treated and charge self-consistency was achieved with net charges 
of ca, +1 on both Che para carbon and nitrogen. These cases closely 
correspond to simple valence bond pictures. The ORTHO case basically 
represents the average of the valence bond structures where first one 
ortho carbon has a net charge of +1 and then the other ortho carbon has 
a net charge of +1. Table 1 lists all the eigenvalues for the four 
forms of nitrobenzene. In each case, the first 23 orbitols arc filled. 
Comparison of the MO energies of the NITRO, ORTHO, and PARA entities 
with the corresponding MO energy for the NORMAL form shows that the 
inequalities in Eqs. (3.8a,b) are obeyed without exception, eigenvalue 
by eigenvalue. That is, in all cases 

I ^i, NORMAL I - I ^i, NITRO I’ I ^i, ORTHO I » I ^i, PARA I' 

where 1^^ pi is the magnitude of the i^th eigenvalue (MO) for the F 
form of nitrobenzene. This inequality can be systematically extended 
when the number of charge centers Increases (all other factors such as 
charge magnitudes being equal). Thus, without exception 

I ^i, NITRO I - I ^i, ORTHO I’ I^IjPARaI’ 

The total energy of the various nitrobenzene forms are tabulated in 
Table 2. Comparison shows the massive energy debt incurred by Imposing 
charge separation. All such structures (NITRO, ORTHO, PAIbV) are much 
less stable than normal nitrobenzene. Figure 1 presents a summary of 
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results. It is interesting that shifts ol electronic density between 
ring and the NOg group are much smaller than shifts within the NO 2 or 
benzene. The distribution of electronic density for both ORTHO and PARA 
structures are in agreement with simple valence bond predictions (though 
the magnitudes are different from simple valence bond expectations) . 

The major difference is that even with net charges of ca, +1 on both 
the nitrogen and ring, the net charges on the oxygens are in the neigh- 
borhood of -h and not -1; the latter is the expected value from simple 
valence bond. The modification to lEH permits us to both mimic charge 
distributions and extract energy changes resulting from such distribu- 
tions. The work also suggests that a hybrid MO-Valence Bond scheme is 
possible whereby the elemeiits of the Hamiltonian matrix in this hybrid 
scheme are defined in the following fashion; 

^F(a),F(b) “ ■'’'*'F(a)^'^F(b)‘^’^ ^ ^F(a) “^FCa) ,F(b) ‘'®f^“'^F(a) ,F(b) ^ f^F(a) 

^F(b)^S’(a),F(b)®F(a),F(b)’ 

if the molecular wavefunctlons for the structures NORMAL, NITRO, ORTHO, 
PARA take the form 


4 ' = 
F 


23 

n ip 

1=1 


i,F- 


Here, F(b) ^ molecular interaction configuration constant, 

EpCa) and are the respective energies of the molecular entities 

F(a), F(b). 
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The Krouecker delta 

^'(a),F(b) ■ 1- 

- 0, F(a) 4 F(b). 

The overlap integral between the electronic density distributions repre- 
sented by F(a) and F(b) For the nitrobenzene cases is given by 


S(a),F(b) “ 
When F(a) « F(b), 


23 

It 

i“l 


•^'*‘'±,F(a)'^i,F(b) 



(4.2) 


^(a),F(a) " %(a)* 

On the other band, when F(a) 4 P(b), 

^’‘F(a),F(b) “ ^S(a),P(b)f^F(a) ’^F(b) ^®F(a) ,F(b) ‘ 

Taking linear combinations of the trial solutions leads to 

the requirement that 

|H - ES| = 0. (4.3) 

The elements of H, F(b)' '^^’^^^fiaration interaction terms 

in the usual sense where electronic correlation is obtained through 
mixing excited state with the ground state approximate wavef unctions . 
Here, structures which differ in charge distribution but are derived 
from the same set of atomic orbitals are allowed to mix. This roughly 
corresponds to a simple valence bond scheme. The full mixing of the 
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four density distributions, as depicted in Fig. 1, was not calculated. 
Instead the simpler cases of (NORMAL, NITRO) , (NORMAL, ORTHO), and 
(NORMAL, PARA) have been treated with “ 1* I’he Hamiltonian 

matrix was transformed using a basis set which dlagonalijjed the overlap 
matrix; formally, the same procedure is used with lEH, except that here 
the elements of the overlap matrix are not between two atomic centers, 
but between two molecular electronic density distributions. For our 
problem. 


11 - ES USn'^u’'’HUS“^y’’ 
<v 'vvD 'V Woin ^ 


El , 

'Xf 


(4.4) 


where 1 is the order 2 identity matrix and U is an order 2 unitary 
matrix 



(4.5) 


which obeys the relationship 


U'U 

f\, f\, 


uu' 

•vx, 


1 . 


(4.6) 


The symbol is the diagonalized overlap matrix. It is found that 

the greater the number of centers of charge separation, the smaller the 
overlap v^ith the normal nitrobenzene configuration and correspondingly 
smaller mixing with the ground state (normal nitrobenzene). In fact, 
mixing ORTHO with NORMAL actually raises the final energy slightly (in 
the seventh decimal place) with respect to the normal nitrobenzene ground 
state energy. The results of tUis exercise fire presented in Table 3. 

It is noted, referring to Table 3, that the NITRO form is the most 
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important variant in that it mixes n»re with the NORMAL entity than 
either ORTHO or PARA. The lEH picture for the isolated nitrobenzene 
molecule, at least, appears to suffice. 

5. Conclusion 

Several modifications which permit iutarnuclear interactions as well as 
electronic interactions to be accounted for explicitly have been intro- 
duced into conventional lEH. The electronic energies obey a variational 
principle as a result of the electronic Interaction modification. The 
modified lEH has been applied to a simple, naive case where it was 
desired to impose charge separation and follow the effects. As a con- 
sequence, since certain gross features of valence bond theory could then 
be mimicked, a hybrid MO-Valence bond scheme was devised which employs 
total molecular electronic density distributions In a fashion similar to 
the way lEH incorporates atomic orbitals. This approach should be a 
useful tool for the organic chemist, since it is capable of quantitative 
evaluation of valence bond-type concepts, which still prove to be 
immensely valuable on a qualitative level. On the other hand, the mod- 
ified lEH is n powerful tool for dealing with charge separation and 
for determining the effect of charge separation in reactions. This 
latter application will be the objective of the accompanying paper. 
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Table 1. The one-electron energies® for four variant total molecular 


electronic density distributions of nitrobenzene 


MO 

NORMAL*’ 

NITRO® 

ORTHO^* 

PARA® 


1 

-37.17 

-31.48 

-30.32 

-30.95 


X 


2 

-30.87 

-29.29 

-26.45 

-27.39 


3 

-30.34 

-26.45 

-25.70 

-26.08 


4 

-26.83 

-25.66 

-22.62 

-24.39 

F 

5 

-26.23 

-25.43 

-22.57 

-22.74 

I 

L 

6 

-22.55 

-20.11 

-18.16 

-18.83 

L 

E 

7 

-20.55 

-20.06 

-17.40 

-18.45 

D 

8 

-18.88 

-17.28 

-15.79 

-16.13 

M 

0 

9 

-18.70 

-16.88 

-14.83 

-15.34 

L 

E 

10 

-17.32 

-14.94 

-12.93 

-13.63 

C 

u 

11 

-16.12 

-14.78 

-12.51 

-13.60 

L 

A 

12 

-15.55 

-14.58 

-12.35 

>13.51 

R 

13 

-15.45 

-14.12 

-11.70 

-12.76 

0 

R 

14 

-15.10 

-13.11 

-11.57 

-11.70 

B 

I 

15 

-14.84 

-12.98 

-11.37 

-11.48 

T 

A 

16 

-14.44 

-12.48 

-10.53 

-11.36 

L 

S 

17 

-13.07 

-12.44 

-10.19 

-10.70 


18 

-13.06 

-11.42 

-10.03 

-10.64 


19 

-12.87 

-11.41 

-9.88 

-10.16 


20 

-12.74 

-10.53 

-9.70 

-10.15 


21 

-12.15 

-9.02 

-8.44 

-8.72 


22 

-11.78 

-8.82 

-8.19 

-8.52 
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Table 1 - Continiwd 


MO 

normal'^ 

NITRO*^ 

ORTHO‘S 

PARA® 



24 

-9.72 

-7.66 

-6.78 

-7.02 


k 

25 

-7.16 

-6.95 

-5.39 

-6.29 



26 

-6.47 

-5.70 

-5.07 

-4.86 

U 






N 

27 

-1.98 

-1.89 

-1.4B 

-1.62 

F 






I 

28 

+10.07 

+9.00 

+8.18 

+8.47 

L 






L 

29 

+10.40 

+9.35 

+8.51 

+8.98 

E 






D 

30 

+15.08 

+13.10 

+12.27 

+12.38 








M 

31 

+16.79 

+15.10 

+13.72 

+13.92 

0 






L 

32 

+20.66 

+20.18 

+17.91 

+18 . 70 

E 






C 

33 

+21.60 

+21.18 

+18.41 

+20.10 

u 






h 

34 

+26.04 

+24.40 

+20.38 

+21.82 

A 






R 

35 

+26.07 

+24 . 63 

+22.98 

+22.02 








0 

36 

+33.44 

+28.75 

+ 26 . no 

+27.58 

R 






B 

37 

+45.01 

+43.03 

+37.27 

+38.05 

I 






T 

38 

+45.38 

+43.71 

+38.07 

+40.50 

A 






L 

39 

+62.15 

+60.17 

+52.45 

+53.25 

C 

1.. 

> 

40 

+62.57 

+61.11 

+53.22 

+57.63 



41 

+87.83 

+85.62 

+74.85 

+78.47 

; 



^In eV. 

^Resulting electronic density distribution using conventional lEH 
(see Fig. 1) . 

*^Ca. +1 net charge maintained on the nitrogen (see Fig. 1). 

^Ca. +1 net charge on nitrogen and ca. +h net charge on both ortho car- 
bons (see Fig. 1). 

*^Ca. +1 net charge on nit rogen and ca. +1 charge on para carbon (see 
Fig. 1). 
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Table 2. Energies of four variant total molecular electronic density 
distributions of nitrobenzene 


a 

Type 

Total energy 
of the 

valence electrons 
<eV) 

Total effective 
Internuclear 
energies 
<eV) 

Total 

effective 

energy 

(eV) 

NORMAL 

-856.8 

+727.3 

-129.5 

NITRO 

-764.1 

+739.9 

-24.2 

ORTHO 

-682.9 

+741.3 

+58.4 

PARA 

-711.3 

+740.7 

+29.4 


®See Fig. 1, 
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Table 3. Overlap, coefficients, and energies for the interaction of 
several variant molecular electronic density distributions with NORMAL® 
nitrobenzene 


(Type, NORMAL) 

Overlap^ 

^SoRMAL’ SyPE^'^ 

Ground-state 
energy (eV)'* 

(NITRO, NORMAL) 

4.72x10“^ 

(1,0009,-0.0237) 

-129.5775 

(ORTHO, NORMAL) 

6.31x10“^ 

(1.0000,-0.0003) 

-129.5028 

(PARA, NORMAL) 

3.26X10“^ 

(1.0001,-0.0017) 

-129,5281 


®The electronic density of NORMAL nitrobenzene was obtained using con- 
ventional lEH (see Fig. 1) . 

^Overlap has been only approximately calculated. The overlap 

23 r . 12 

•^'^'tYPe'^NORMAL^'^ “ ^f^[‘^'*’l,TYPE'*’i, NORMAL '*’i,I 

filled molecular orbital in an electronic configuration described by 

'I'l- 

®The state “ *^NORMAl\oRMAL Sype'^TYPE* 

j 

The ground-state energy of NORMAL nitrobenzene “ -129.5028 eV. 
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FIGURE CAPTION 


Fig. 1, The net charges and total overlap populations between neighbor- 
ing atoms for four variants of the electronic density distribution for 
nitrobenzene. The structure labelled NORMAL was calculated with lEH. 
NITRO, ORTHO, and PARA used the electronic interaction modification to 
achieve a particular net charge on selected atoms, e.g. , the nitrogen, 
ortho, and para carbons. The NITRO structure achieved charge self- 
consistency with the net charge shown on the nitrogen (0.97). ORTHO 
achieved charge self-consistency with the net charges shown on the 
nitrogen (0.91) and the ortho carbons (0.54). PARA achieved charge 
self-consistency with the net charges shown on the nitrogen (0.94) and 
the para carbon (0.81). The number in parentheses in the middle of a 
line on the left side of a structure is the total overlap population 
between the two atoms connected by that line. All other numbers (center 
and right) are the net charges on the appropriate atoms 



